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We investigate the ground state properties of an 5=1/2 distorted diamond chain described 
by the Hamiltonian H = Ji J2i{(^st-i • S M + S 3e ■ S 3e +i) + J2SU-2 ■ S M -\ + -h (5 3f _ 2 ■ S M + 
S31 ■ S-ii+i) —HSg} (Ji, J2, J3 >0), which models well a trimerized 5 = 1/2 spin chain system 
Cu3Cl6(H20)2-2HgC4S02. Using an exact diagonalization method by means of the Lanczos 
technique, we determine the ground state phase diagram in the H = case, composed of the 
dimerized, spin fluid, and ferrimagnetic phases. Performing a degenerate perturbation calcu- 
lation, we analyze the phase boundary line between the latter two phases in the J2 , J3 <C Ji 
limit, the result of which is in good agreement with the numerical result. We calculate, by 
the use of the density matrix renormalization group method, the ground state magnetization 
curve for the case (a) where Ji = 1.0, J2 = 0.8, and J3 = 0.5, and the case (b) where Ji = 1.0, 
J2 =0.8, and J3 =0.3. We find that in the case (b) the 2/3-plateau appears in addition to the 
1/3-plateau which also appears in the case (a). The translational symmetry of the Hamiltonian 
Ti. is spontaneously broken in the 2/3-plateau state as well as in the dimerized state. 
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§1. Introduction 

In the past years a great deal of work has been devoted 
to the study of quantum spin systems with competing 
interactions. This has stemmed largely from their be- 
ing systems which include two phenomena of great inter- 
est, frustration and quantum fluctuation. Recently, mag- 
netic properties of a trimerized 5=1/2 spin chain sys- 
tem Cu3Cl6(H20)2-2H8C4S02 have been experimentally 
studied by Ishii, Tanaka, Hori, Uekusa, Ohashi, Tatani, 
Narumi, and Kindo. 1 ) (See below for the discussion on 
the competition of the exchange interactions in this sys- 
tem.) They have measured the temperature dependence 
of the magnetic susceptibility and also the magnetiza- 
tion curve at sufficiently low temperatures. Their re- 
sults clearly demonstrate that the ground state of this 
system is nonmagnetic and that there exists the energy 
gap with the magnitude of about 3.9 T or 5.2 K between 
the ground state and a first exited magnetic state. Fur- 
thermore, the result for the magnetization measurement 
shows that the so-called 1/3-plateau in the magnetiza- 
tion curve may start from about 58 T. 

As has been discussed by Ishii et al., l > the magnetic 
properties of Cu3Cl6(H 2 0)2-2H 8 C4S02 can be well de- 
scribed by the following Hamiltonian TL, which is a sum 
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of the exchange interaction term 7i GX and the Zeeman 
term 7ii: 

H=H CX + H Z , (1.1a) 

AT/3 

T~(-cx = Ji /J (631-1 • S 3 e + 5 3 £ • S 3 i + i) 

1=1 

N/3 

t=\ 

N/3 

+ J3 ^](5'3£-2 ' S3I + S31 ■ SW+2) ) (1-lb) 
l=\ 
N 

h z = -hJ2s!, (Lie) 
1=1 

where Si is the spin operator with the magnitude 5=1/2 
for the Cu 2+ ion located at the £th site, and periodic 
boundary conditions (Sn+i = 5i,5jv+2 = £2) ar e im- 
posed. In the following discussions, the total num- 
ber N of spins is assumed to be a multiple of six. 
From the crystal structure and the lattice parameters 
of Cu 3 Cl 6 (H 2 0)2-2H 8 C4S02, Ishii et al.^ have inferred 
that three exchange constants, J\, J2, and J3 are all an- 
tiferromagnetic (positive). Thus, the three interactions 
compete with each other. They have also inferred that 
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Ji is larger than J2 and J3. This means that the spins 
>Sa«_i, S 3 i, and S 3 e+i form a trimer; J\ is the intra- 
trimer exchange constant, and J2 and J3 are the inter- 
trimer ones. Schematical representations of the Hamil- 
tonian H cx are given in Figs. 1(a) and 1(b). Referring to 
Fig. 1(a), we may call this model the 'distorted diamond 
chain model'. 




Fig. 1. (a) and (b) Two ways of representing schematically the 
Hamiltonian Hex- (c) and (d) The two-fold degenerate dimerized 
(D) states. In (a) -(d), the open circles stand for the S = 1/2 
spins, and the solid, wavy, and dotted lines correspond to the 
Jl, J2, and J3 interactions, respectively. Two spins surrounded 
by an ellipse in (c) and (d) form a singlet dimer 



In the present paper, we explore the ground-state 
properties of the distorted diamond chain, performing 
mainly numerical calculations based on the Lanczos 
method as well as the density matrix renormalization 
group (DMRG) method. 2 ' 3 ^ We focus our attention upon 
the ground state phase diagram in the case of zero ex- 
ternal magnetic field (H = 0) and the ground state mag- 
netization curve. We treat the case where J\ > J3 and 
J2 > 0, noting that the J% < J3 case is equivalent to the 
Ji > J3 case by interchanging the role of J\ and J3 . 

Let us now discuss three special cases where H = is 
assumed. First, when J2 = 0, we can readily show by the 
use of the Lieb-Mattis theorem 4 ) that the ground state 
is a ferrimagnetic (FRI) state where the magnitude .Stot 
of the total spin S to t=J2e=i & O^tot = 'Stot (Stot + 1)) is 
given by Stot = N/6. We note that this value is one third 
of that of the maximum (ferromagnetic) value, N/2. Sec- 
ond, when J3 = 0, the system is reduced to the J\-J\-Ji 
trimerized 5 = 1/2 antiferromagnetic chain; in particular, 



if Ji = J 2 , the 5 = 1/2 uniform antiferromagnetic chain is 
obtained. Thus, the ground state is the spin fluid (SF) 
state with Stot = 0, which has no excitation energy gap. 5 ) 

Third, the case where J\ = J3 has already been inves- 
tigated by Takano, Kubo, and Sakamoto. 6 ) They have 
shown that the ground state is the FRI state with S tot = 
iV/6, the tetramer-dimer (TD) state with 5 tot = 0, or the 
dimer-monomer (DM) state with Stot^O, depending on 
whether J 2 /Ji < 0.909, 0.909 <J 2 /J\ < 2, or 2 < J 2 / J x . 
In the TD state, quadruplets £3^-3, £3^-2, £3^-1 , an d 
S 3 i (or S 3e , S 3e+1 , £3^-2, and £3*4-3) of spins form singlet 
tetramers, and pairs S 3 £+i and £3^+2 (£3^4-4 and £3^+5) 
of spins, which we call J 2 -pairs, singlet dimers. Thus, 
the TD state is two-fold degenerate. In the DM state, 
on the other hand, the J2-pairs form singlet dimers, and 
the remaining monomer spins S 3 g are free, which leads 
to 2 JV / 3 -fold de generacy. These are schematically illus- 
trated in Fig. 1 in ref. 6. The explicit expressions of the 
eigenfunctions and some physical quantities for both the 
TD and the DM state are given in Appendix. As has 
been shown by Takano et aZ., 6 ) the commutation rela- 
tion \H CX , (Su+i + S , 3£ +2 ) ] = holds in the Ji = J 3 
case. Therefore, each J2-pair forms a singlet or a triplet. 
It is interesting to note that in the DM and FRI states, 
all the J2-pairs form the singlets and the triplets, re- 
spectively, while in the TD state the singlet and triplet 
J2-pairs are arranged alternatively. 

In a recent paper, 7 ) Okamoto, Tonegawa, Takahashi, 
and Kaburagi have studied the ground state of the dis- 
torted diamond chain in the case of H — 0, employing an- 
alytical methods based mainly on the bosonization tech- 
nique and doing physical considerations. According to 
them, a tetramer in Takano et aZ.'s TD state is very pe- 
culiar to the Ji = J3 case, and when J\> J 3 , the tetramer 
consisting of the spins, £3^-3, £3^-2, £3^-1, and S 3 i is de- 
composed into two dimers, one of which consists of £3^-3 
and S 3 t-2, and the other of £3^-2 and £3^—1. Thus, 
the TD state is a special case of the two-fold degener- 
ate dimerized (D) state with £ tot = 0, shown in Figs. 1(c) 
and 1(d). The DM state is also peculiar to the J\ = J 3 
case. When J\> J 3 , the monomer is no longer free and 
interacts with neighboring monomers through the dimer 
between them. Thus, the DM state is a special case of 
the SF state with St o t = 0. 

Summarizing these, Okamoto et aL 7 ) have concluded 
that the ground state phase diagram of the distorted di- 
amond chain is composed of the D and SF phases plus 
the FRI phase which is stable at least when J2 is suffi- 
ciently smaller than Ji and J3. They have also shown 
that the transition between the D and SF phases is of 
the Berezinskii-Kostclitz-Thouless (BKT) type, as in the 
case of the S = 1/2 antiferromagnetic chain with uni- 
form nearest and next-nearest neighbor interactions, 8 ) 
and have presented the ground state phase diagram de- 
termined numerically. In the region where the D state 
is the ground state, there exists a finite gap between the 
two-fold degenerate ground state and a first exited state. 
As can be seen from Figs. 1(c) and 1(d), the period of 
the translational symmetry of the D state is six, which is 
twice as large as that of the Hamiltonian 7i ex . Thus, the 
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spontaneous symmetry breaking occurs in the D ground 
state. This is consistent with the necessary condition 
for the appearance of the plateau in the ground state 
magnetization curve, 

n(S — m) = integer , (1.2) 

obtained by Oshikawa, Yamanaka, and Affleck, 9 ) Here, n 
and m are, respectively, the period of translational sym- 
metry of the plateau state and the average magnetization 
per site in the plateau, and S is the magnitude of spins 
in the system. (Note that n = 6, S = 1/2 and m = for 
the present ground D state in the H = case.) 

In the present paper, we discuss the derivation of the 
ground state phase diagram in more detail. Further- 
more, we discuss the ground state magnetization curve, 
as mentioned before. 

§2. Ground State Phase Diagram in the H = 
Case 

Let us denote, respectively, by (Stot', N) and 
E^(Stot', N) the lowest and second-lowest energy eigen- 
values for a given set of 5 to t and N of the Hamiltonian 
H cx . We have calculated numerically these eigenvalues 
for finite size systems with N = 6, 12, 18, and 24, employ- 
ing the computer program package KOBEPACK/S 10 ) 
coded by one of the present authors (M. K.) by means 
of the Lanczos technique. Our calculation shows that, 
depending on the values of Ji, J2, and J3, only 
2?(°)(0; N) or E(°\N/6; N) becomes minimum among 
E^(Stof,Ny S for all values of S tot (= 0, 1, N/2). 
This is consistent with the fact that the ground state 
phase diagram is composed of the D, SF, and FRI 
phases, 7 ^ as discussed in § 1. 

The ground state phase is the FRI phase when 
£(°)(AT/6; oo) < E (a \0; oo), and it is the D or the SF 
phase when E^(N/6; oo) > £(°)(0; oo). We estimate 
the critical value J™ 1-0 of Jj (i = 1 2, or 3) between the 
former phase and one of the latter phases for a given set 
of the other two J's in the following way. We first calcu- 
late numerically the value J™ 1 " (A) which satisfies 

E^(N/6; N) = £(°)(0; N) . (2.1) 

Then, we estimate J™~° by fitting J F c RI ~°(A)'s with 
12, 18, and 24 to a quadratic function of l/N 2 , that is, 

J™"°(A0 = J™-° + £ + ^ , (2.2) 

where a and b are numerical constants. 

On the other hand, in order to estimate the critical 
value Jj D c ~ between the D and SF phases, we employ 
the method of level spectroscopy, 8 ' u ) which has been 
successfully applied to estimate numerically the BKT 
critical points in many cases. The procedure of the 
method is as follows. (For the physical interpretation 
of this procedure the reader is referred, for example, to 
ref. 7.) First, we introduce the singlet-singlet energy gap 
A SS (A) and the singlet-triplet energy gap A SS (A) for the 
finite- N system defined by 

A ss (N) = (0; N) - (0; N) , (2.3a) 



A rt (N) = £ (0) (1; N) - £(°>(0; N) , (2.3b) 

where E^ (0; N) is nothing but the ground state energy. 
Then, we calculate the value </j D c ~ SF (A) satisfying 

A SS (A) = A st (AO • (2.4) 

Finally, we extrapolate the results jf~ SF (N)'s to N —* oo 
by the use of a polynomial of 1/A 2 to estimate Jf~ SF - 
In the practical calculations we have performed this ex- 
trapolation by using the values of J°T SF (A)'s with 12, 
18, and 24 and the formula, 

J° C - SF (N) = J°- SF + ^ + ^ (2.5) 

with numerical constants a' and b' . 

We can obtain the ground state phase diagram on 
a exchange constant parameter plane by plotting the 
estimated critical points and by connecting them as 
smoothly as possible. The results are depicted in 
Figs. 2(a) and 2(b), which show, respectively, the phase 
diagram on the J3 versus J2 plane with J\ fixed at J\ = 1 
and that on the J3 versus J\ plane with J\ = 1 ; note that 
the former is a rcplot of Fig. 9 in ref. 7. In these phase di- 
agrams the case where Ji < J3 is included to make them 
more complete. 

The phase transition between the FRI phase and the 
D or the SF phase is, of course, of first order. Further- 
more, when Ji = J3, the transition between the D and 
SF phases, or cquivalcntly, that between Takano et al.'s 
TD and DM phases, at J2/J1 = 2 is also of first order. 
As is shown in Appendix, the ground state energy e g (N) 
per site for 0.909 < J2/J1 < 2 and that for the 2<J 2 /J\ 
case are, respectively, given by s g (N) = Etd(N)/N = 
-(4Ji + J 2 )/12 and e s (N) = E UM {N)/N = - J 2 /4, and 
therefore, plotting de g {N)/dJ\ as a function of Ji, we 
have a discontinuity at J\ = Ji/2. 

We assert, for the following reasons, that the first or- 
der phase transition is peculiar only to the above two 
cases in the present system. Let us first give a physical 
consideration. It is well known that the ground state of 
the S—l/2 chain with uniform, antifcrromagnetic near- 
est and next-nearest neighbor interactions is a two-fold 
degenerate dimerized state which is essentially equiva- 
lent to those shown by Figs. 1(c) and 1(d), as far as the 
latter interaction is sufficiently larger than the former 
one. 8, 12, 13 -* If a bond alternation is introduced in the 
nearest neighbor interaction as a perturbation, the de- 
generacy of the ground state is lifted, and one of the 
dimerized state or the other becomes the ground state 
depending upon the value of a bond alternation param- 
eter. Thus, sweeping this parameter, we have the first 
order transition in the ground state at the uniform point. 
In fact, the first order transition has been observed, for 
example, in the ground state of a generalized S = 1/2 
ladder with additional exchange interactions on diago- 
nal bonds, 14 ) which has bond alternating nearest neigh- 
bor interactions when represented in a one-dimensional 
chain scheme like Fig. 1(b). In the present distorted di- 
amond chain, on the other hand, we have no bond al- 
ternation, but have a bond trimerization in the nearest 
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(and also in the J\ > J3 case) except for the FRI-D tran- 
sition. 
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Fig. 2. Ground state phase diagrams (a) on the J3 versus J2 
plane with J\ = 1 and (b) on the J3 versus J\ plane with J\ = 1. 
The phase diagrams are composed of the dimerized (D), spin 
fluid (SF), and ferrimagnetic (FRI) phases. 



neighbor interaction [see Fig. 1(d)] which never lifts a 
two-fold degeneracy in the D ground state. Thus, no 
first order transition occurs in the ground state of the 
present chain, except for the two peculiar cases stated 
above. 

In order to supplement the above reason discussed 
physically, we have calculated numerically the J\- 
dependence of de E {N)/dJi with iV = 6, 12, 18, and 24 
for the J2 = 1 and J 3 / J x = 0.999 case, for which the D-SF 
critical point is given by Jic 0.5002 according to the 
phase diagram shown in Fig. 2(b). The result is plotted 
in Fig. 3, together with that for the J2 = 1 and J3/J1 = 1 
case mentioned above. This figure demonstrates clearly 
that de e (N)/dJ 1 for the J 2 = 1 and J3/J1 =0.999 case 
have no anomaly around J\ = 0.5. This result suggests 
that we have no first order transition in the J± < J3 case 



Fig. 3. Plot versus J\ of de g (N)/dJi with N = 6 (closed circles), 
12 (open circles), 18 (pluses), and 24 (crosses) for J2 = 1 and 
J3/J1 = 0.999, together with that (solid lines) for J2 = 1 and 
Ja/Ji = l. 



§3. Ground State Magnetization Curve 

We have calculated the ground state magnetization 
curve of the present system, using the DMRG method 
proposed originally by White. 2 -' Here, we describe briefly 
the procedure for our DMRG calculation. We employ the 
finite system algorithm 3 - 1 as improved by White, 15 -* which 
reduces substantially the required computational time. 
The maximum number of block states we keep in the 
calculation is m max = 100, which leads to the truncation 
error of the order of 10~ 8 . Comparing the calculated re- 
sults for the cases where m = 70 and m= 100 block states 
are kept, we may conclude that the results for m=100 
are within an accuracy of roughly 10~ 4 ~ 10~ 5 . Thus, 
we adopt these results as our final results without carry- 
ing out any m-extrapolation. We assume open boundary 
conditions for the Hamiltonian H cx , subtracting the two 
terms, J\ Sn ■ Sn+i(= JiSn ■ Si) and J3 Sn ■ 5jv + 2(= 
J3SN ■ S2), from the expression of eq. (1.1b). This is be- 
cause, as is well known, open boundary conditions make 
the DMRG method work more effectively. It should be 
noted that these boundary conditions lead to an open 
chain which has no inversion symmetry with respect to 
its center. 

For finite size systems with up to iV = 96 spins, we have 
calculated the lowest energy eig envalue E^(M; N) of 
7i cx within the subspace determined by the value M of 
the z-component of S^ot- Once the values of E^(M; N) 
for nonnegative M's (M — 0, 1, • ■ ■, N/2) are known, the 
ground state magnetization curve can be obtained by 
plotting as a function of H the average magnetization 
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m(N)(=M/N) per site determined from 



m{N) 







with M = l, 2, 



(when AE^ 0) (M, M - q; N) > H) , 

max^M I AE m (M,M-q; N)<H 1 , 

— 1 (otherwise) , 

(3.1) 

M, where 



JV 



, AT/2, and g=l, 2, 



E^(M-q; N) 



AE^(M,M-q; N) 

y 

(3.2) 

It is noted that m(N) = when < H < H (N) , where 
H (N)=mm{AE^ 0) (q,0; N)}, and the saturation field 
H a (N) is given by # s (iV) = max{A^ 0) (f , f - 1; JV)}. 
The resulting magnetization curve is a stepwisely in- 
creasing function of H. In particular, when the 
ground state in the case of H = Q is the Stot = 
state and the competition among the three interac- 
tions is not too strong, the magnetization curve has 
N/2 steps, starting from m(N) — and ending at 
m(N) = 1/2; at each step m{N) increases by 1/7V. 16 ) 
Then, H {N) = A£[°\l,0; N) = A at {N), and H a (N) = 
AE[ \M,M-1;N). Following Bonner and Fisher's pi- 
oneering work, 17 ) we may obtain, except for plateau re- 
gions, a satisfactorily good approximation to the magne- 
tization curve in the thermodynamic (TV— >oo) limit by 
drawing a smooth curve through the midpoints of the 
steps in the finite size results. 



Table I. Numerical values of H (N)/H a (N), H 1/3 ^^/H^N), 
H 1 /3,hW/-HsW, H 2/3il (N)/H s (N), and H 2/3>h (N)/H s (N) 
estimated in the cases (a) and (b). Note that the figures in the 
parentheses for the N-^oo values show errors in the last digit. 



(a) the case 


(a) where Ji 


= 1.0, J 2 = 0.S 


S, and J3 = 0.5. 






Af = 48 


7V = 72 


Af = 96 




H (JV) 
H S (N) 
Hl/3.lW 
H B {N) 

HJN) 


0.024525 


0.016736 


0.012795 


0.0033 (5) 


0.160260 
0.702154 


0.160875 
0.699845 


0.161043 
0.699090 


0.1612 (1) 
0.6982 (1) 


(b) the case 


(b) where Ji 


= 1.0, J 2 = 0. 


3, and J3 = 0.3. 






Af = 48 


iV = 72 


JV = 96 


N~>oo 


H (N) 

H B (N) 
H 1/3A (N) 

HJN) 
Hl/3,hW 

H S (N) 


0.031482 
0.312021 


0.021075 
0.314803 


0.015841 
0.315499 


0.0 
0.3162(1) 


0.842373 


0.841792 


0.841618 


0.8414(1) 


%/3,l(*) 

H S {N) 
H 2/3 ,h(JV) 


0.902840 


0.905461 


0.906667 


0.9084 (1) 


0.944329 


0.942955 


0.942391 


0.9416 (1) 



Figure 4 shows the magnetization curves calculated 
for N — 96; Figs. 4(a) and 4(b) are, respectively, for the 
case (a) where J\ = 1.0, J2 = 0.8, and J3 = 0.5, and 
for the case (b) where J\ = 1.0, Ji = 0.8, and J3 = 0.3. 
In each case, we have a 1/3-plateau, and moreover, in 
the case (b), we have a 2/3-plateau in addition to this. 
We denote the highest and lowest values of H giving the 
p-plateau (p=l/3, 2/3) by H p ^(N) and H p ^(N), respec- 




0.5 1.0 

H/H s 




0.5 1.0 

H/H s 



Fig. 4. Ground state magnetization curves calculated for TV = 96 
(a) in the case (a) where Ji = 1.0, J2 = 0.8, and J3 = 0.5 and (b) 
in the case (b) where Ji = 1.0, J2 = 0.8, and J3 = 0.3. 



tively. We have calculated these as well as Hq(N) and 
H S (N) also for TV = 72 and 48 in the cases (a) and (b), 
and have extrapolated the results to N — > oo to estimate 
the values in the thermodynamic limit. The results are 
tabulated in Table I. We see from this Table that, except 
for Ho(N), the finite size results for iV = 96 give good ap- 
proximations to the results in the thermodynamic limit. 

The 1 /3-plateau state is essentially the same as the fer- 
rimagnetic state which becomes the ground state when 
J 3 >J 3 ™" - 0.813 in the J 1 = 1.0, J 2 = 0.8, and H = 
case. On the other hand, we see from eq. (1.2) that 
the periodicity n of the 2/3-plateau state should be a 
multiple of six, since ,5=1/2 and M = l/3. In order to 
check this fact, we have calculated, by means of the 
DMRG method, the ^-dependence of the expectation 
value m(£; N) of S| for the 2/3-plateau state in the 
case (b). The results for TV = 96 depicted in Fig. 5 clearly 
demonstrates that n = 6 for this state, which means that 
the translational symmetry is spontaneously broken also 
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Fig. 5. Plot versus I of m(l; N) with N = 96 for the 2/3-plateau 
state in the case (b) where Ji = 1.0, J2 = 0.8, and J3 = 0.3. Solid 
lines are guides to the eye. 



in this state as in the D state. Thus, the mechanism 
for the appearance of the 2/3-plateau in the present 
case is considered to be the same as that, which has 
been clarified by Totsuka, 18 ' for the appearance of the 
1/2-plateau in the ground state magnetization curve of 
the 5=1/2 antiferromagetic chain with bond alternat- 
ing nearest and uniform next-nearest neighbor interac- 
tions. 19 ^ Based on this mechanism, we can determine, by 
using the method of level spectroscopy, 8, n ' the region, 
where the 2/3-plateau appears, of Ji for a given set of 
the other two J'. Leaving the details of this investigation 
for a forthcoming paper, we only mention here that, ac- 
cording our preliminary result, the 2/3-plateau appears 
when 0.169 < J 3 < 0.375 in the J l = 1.0 and J 2 = 0.8 case. 

§4. Concluding Remarks and Discussion 

We have explored the ground state properties of 
an .5=1/2 distorted diamond chain described by the 
Hamiltonian Ti [see eqs. (l.la)-(l.lc)], which mod- 
els well a trimerized S=l/2 spin chain system 
Cu3Cl6(H20)2-2HgC4S02. Using an exact diagonaliza- 
tion method by means of the Lanczos technique, we have 
determined the ground state phase diagram in the H = 
case [Figs. 2(a) and 2(b)], composed of the D, SF, and 
FRI phases. As has been mentioned in §1, the ground 
state of Cu3Cl6(H 2 0)2-2H 8 C4S02 is nonmagnetic with a 
relatively small but finite energy gap. Furthermore, the 
three exchange constants, J±, J2, and J3, in this system 
satisfy J\ > J2, J3 > 0. From these facts together with 
the obtained phase diagram shown in Figs. 2(a), we may 
anticipate that in Cu3Cl6(H20)2-2HgC4S02 the values 
of J 2 / Ji and J3 / Ji are in the region of 0.7 < J 2 / Ji < 0.9 
and 0.45 < J3/ Ji < 0.65. We have also calculated, by the 
use of the density matrix renormalization group method, 
the ground state magnetization curve for the case (a) 
where J x = 1.0, J 2 = 0.8, and J 3 = 0.5 [Fig.4(a)], and the 
case (b) where Ji = 1.0, J 2 = 0.8, and J 3 = 0.3 [Fig. 4(b)]. 
We have found that in the case (b) the 2/3-plateau ap- 
pears in addition to the 1/3-plateau which also appears 



in the case (a), and have clarified that the translational 
symmetry of the Hamiltonian 7i cx is spontaneously bro- 
ken in the 2 /3-plateau state. 

Figure 2(a) demonstrates that the boundary line be- 
tween FRI and SF phases in the limit of < J2 <C Ji 
and < J3 <C Ji is approximately given by J 2 = J3 • This 
can be understood analytically by performing a degener- 
ate perturbation calculation around the truncation point 
J2 = J3 = in the following way. The unperturbed Hamil- 
tonian 7io is given by the sum of the Hamiltonian he for 
the £th trimer consisting of S 3 g~\, S 3 e, and S 3 e+i: 



N/3 

TCa = ^2 he . 



he — Sat-i • S 3 e + S^i ■ S 3 e+i ■ 



(4.1a) 
(4.1b) 



where the value of Ji is taken to be Ji = 1.0 as the unit 
of energy. The eigenfunctions, 4>P and 4>j, of the lowest 
energy states of he, which are two-fold degenerate, are 
expressed, by the use of ae and fle representing, respec- 
tively, the Sf = 1/2 and S\ = —1/2 single spin states, 
as 

<fiP = -j={\otU-iOtuPu+x) - 2 



V6 1 



t+i) 



(.(2) 



V6 



+ \/3 3 e-ia 3 ea 3 e+i}} , (4.2a) 

{\PM-lp3e°l3t+l) — 2 |/33£-ia3£/33£+l) 

+ \a 3 e-i(3 3 e(3 3 e+i)} ■ (4.2b) 



Restricting only to these two eigenfunctions, we can rep- 
resent S 3 e-i, S 3 e, and S 3 e+i in terms of the pseudo 
5 = 1/2 operator Te associated with the £th trimer: 



J 3£-l 



-To 



QZ rpZ 



°3£ ~ 3 L l 



'->>, 



s. 



3£+l 



= ~n 



~iZ rpZ 

>3l+\ — ~^ L l i 



(4.3a) 
(4.3b) 
(4.3c) 



where the Tf = 1 /2 and T/ = — 1/2 states correspond to 
4>P and 4>l , respectively. Substituting eqs. (4.3a)-(4.3c) 
into the second (J 2 ) and third (J3) terms in the right- 
hand side of eq. (1.1b) leads to the following effective 
Hamiltonian Ti. c s : 



N/3 



W e ff — - {J 2 - J3) ^2 ^ e ' Tz+l 



(4.4) 



with Ti 



(N/3)+r- 



zT±. We see from this equation that the 



ground state of the pseudo Te spin system is the fer- 
romagnetic or the SF state depending upon whether 
J 2 > J3 or J2 < J3, the former and the latter correspond- 
ing, respectively, to the FRI and SF states in the original 
Se spin system. Thus, it is shown that the J 2 = J3 line 
yields the boundary line between FRI and SF phases in 
the limit of < J 2 < Ji and < J 3 < Ji . 
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Appendix 

Here, we discuss the Ji = J 3 case. 6 ) Let us introduce 
the following states: 

[£,£+1,1 + 2,1 + 3] 

{\a>£ai + if3£ +2 (3i +3 ) + \f3il3i +1 ai +2 ai +3 ) 



+ \ai(3i + ia e+2 (3i +3 ) + \(3 e a e+1 (3 e+2 a e+3 ) 
- 2 \aePe + i/3t + 2ae +3 ) - 2 \/3tae + iat+2Pe+a)} , 

(A-la) 



[£,£+!} = — {|a*ft+i> -|&a/+i>} > 



(Alb) 
(Ale) 



Then, the eigenfunctions <&^ D (A) of the TD state for the 
finite- AT system is given by 



1 



$± B (N) = ^{^l(N)±$W(N)} 



(A2) 



with 
$«(A) 



[1,2] [3, 4, 5,6] ■■■ 

[N — 5, A— 4] [AT— 3, A-2, A- 1, A] 



(A-3a) 



(2), 



, TD (A) 

= [N, 1,2, 3] [4, 5] ••• 

[A- 6, A- 5] [A— 4, A- 3, A- 2, A- 1] . 

(A-3b) 

Note here that (^(JV)|$^,(JV)) = 0. On the other 
hand, the eigenfunction $dm(A) of the DM state for 
the finite-A system is 

$dm(A) = [1, 2] [3] ••• [A- 3, A- 2] [A] . (A-4) 

It is easy to show that the energies of the TD and DM 
states for the finite-A system are given, respectively, by 
£td(A) = -$(4Ji + J 2 ) and E DM (N) = -f J 2 . We 
note that Eum(N) is the total energy of y indepen- 
dent singlet dimers mentioned in § 1, while Etd(N) is 
the total energy of ^ independent singlet tetramers plus 
that of independent singlet dimers. Furthermore, it is 
also straightforward to calculate the spin pair correlation 
functions in the TD and DM states, defined, respectively, 



by 

WTDit,?; N)= ($± D (A)|5, • S e ,\$± D (N)) , (A5a) 

udm(£,£'; N)= ($± M (N)\S e ■ S e ,\$± M (N)) . (A5b) 

When A = 12, 18, • • •, the function uj T r>(£,£'; A) for £< 
£' 20 ) takes nonzero values independently of A only in 
the following cases: 



utu(£,£; N)=-, 



w T d(3^ - 2, 'S£; A) = w TD (3^, U + 2; A) 
1 

~ ~12 ' 

w TD (3^- 1,3^; N)= lo td (3£,3£+ 1; A) 
1 

~ ~12 ' 
1 



(A-6a) 

i 

(A-6b) 



wtd(3£ + l,3£ + 2; A): 
w T d(3^ + 3; A): 



12 



1 

24 



(A-6c) 
(A-6d) 
(A-6e) 



On the other hand, the nonzero values of u>dm(£,£'', N) 
for £<£' 2 °) are 



ujbm{£,£; A)= - , 



wdm(3*+1,3* + 2; A) = 



(A-7a) 
(A-7b) 



which are independent of A(=6, 12, 18, • • •). 
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